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Abstract 

We calculate the tunneling density of states (DOS) of MgB2 for different tunneling directions, by directly solving 
the real-axis, two-band Eliashberg equations (EE) . Then we show that the numeric inversion of the standard single- 
band EE, if applied to the DOS of the two-band superconductor MgB2, may lead to wrong estimates of the strength 
of certain phonon branches (e.g. the Ei g ) in the extracted electron-phonon spectral function a 2 F(uj). The fine 
structures produced by the two-band interaction turn out to be clearly observable only for tunneling along the ab 
planes in high-quality single crystals. The results are compared to recent experimental data. 



Key words: electron-phonon interaction. Eliashberg equations, magnesium diboridc 



There is a growing consensus that superconductivity 
in MgB2 is driven by the electron-phonon interaction 
(EP1) [1] . The idea of multiband superconductivity in 
MgB2 [2,3,4,5] is supported by many recent experi- 
mental results from tunneling [6,7], point contact [8,9] 
and specific heat measurements [10]. These data give 
direct evidence that the superconducting gap has two 
different values: A CT for the two quasi-two-dimensional 
a bands, and A T for the 3D 7r bands [2,3]. According 
to most calculations [2,4,11], the EP1 or, equivalently, 
the Eliashberg spectral function a 2 F(uj) (EF) is domi- 
nated by the optical boron bond-stretching E^g phonon 
branch around 60 - 70 meV. 

At present, the main experimental tool for the de- 
termination of the EPI in superconductors is tunnel 
spectroscopy. Unfortunately, the standard single-band 
procedure to obtain the EF from the first derivative 
of the tunneling current is restricted to a momentum- 
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independent s-wave order parameter and cannot be 
used in anisotropic superconductors like MgB2 . Never- 
theless, there has been a recent attempt to obtain the 
EPI in MgB2 by using this standard approach [12]. 

In the present paper we clarify what information can 
be obtained by using this single-band standard proce- 
dure if one applies it to a two-band superconductor. 

The parameters for the two-band model used in this 
work are based on first-principle electronic structure 
calculations [11]. The four interband and intraband 
electron-phonon spectral functions o?Fij(uj) (where 
i, j — it, a) and the Coulomb pseudopotential matrix 
are the basic input for the two-band Eliashberg 
theory. The aa EPI is dominated by the optical boron 
bond-stretching E2 g phonon mode. For the other 
channels there are also important contributions from 
low-frequency (30-40 meV) and from high-frequency 
phonon modes (~ 90 meV). 

We obtained the theoretical conductance curves 
of MgB2 for different tunneling directions by solv- 
ing the corresponding two-band Eliashberg equations 
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Fig. 1. The calculated tunneling DOS at T— K in the a&-planc 
(a) and along the c-axis (b). The two insets show the fine 
structures of the tunneling DOS due to the electron-phonon 
interaction. 

(EE) [3,13] in the real-axis formulation. The only 
free parameter is the normalization constant \i in 
the Coulomb pseudopotential matrix which is fixed 
in order to reproduce the experimental T c = 39.4 K. 
Figures 1 (a) and (b) show the calculated conduc- 
tances for tunneling along the ao-plane and the c-axis 
directions, respectively. According to Ref. [3], these 
conductances are a weighed sum of the a- and 7r-band 
contributions to the DOS, whose weights ui CT and w n 
are determined by the corresponding Fermi velocities 
in the bands and by the angle of the tunneling cur- 
rent with respect to the ab plane. The values of these 
weights are w a =0.33 and «;^=0.67 for tunneling along 
the ab plane, w CT =0.01 and and u>^=0.99 for tunneling 
along the c axis. It is clear that the contribution of the 
7T band is always dominant even if tunneling is almost 
in the ab plane. The fine structures due to the EPI 
have maximum amplitudes of the order of 0.5% and 
2-3% for measurements along the c-axis and in the 
afe-plane, respectively (see the two insets of Fig. 1). 

In order to test the reliability of the numeric inver- 
sion of the standard single-band Eliashberg equations 
when applied to a two-band superconductor, we used 
the calculated tunneling DOS shown in Fig. 1 as an in- 
put. Of course, in this case, the inverted spectral func- 
tion should correspond to a mixture of the a- and n- 
band contributions. Figures 2 (a) and (b) show the re- 
sults of the inversion of the calculated tunneling DOS 
of Fig.l (a) and (b), respectively, using a standard 
single-band code (solid lines). The effect of the single- 
band EE inversion is estimated by a least-square fit 
of the inverted spectral functions with the weights for 



a 2 F aa + a 2 F a 



and a 



as free parameters. These fits are shown by the dashed 
curves in Fig. 2 and correspond to ot 2 F^ ~ 0.31a 2 F a + 
0.160^ and a 2 F/''* ~ 0.01a 2 F a + 0.99a 2 F n . 

Tunneling along the c-direction practically corre- 
sponds to a single-band case and, as expected, the 
inversion properly reproduces the a 2 F^. In contrast, 
in the case of tunneling along the afe-direction, the a- 
band spectral functions play an essential (and ampli- 
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Fig. 2. The Eliashberg functions for tunneling along the 
arj-plane (a) and along the c-axis (b) obtained from the curves 
of Fig.l by inversion of the single-band Eliashberg equations 
(solid line). The dashed lines represent the least-square fits 
with two free parameters as explained in the text. 

fied) role since about 33% of the cr-band contribution 
in the tunneling DOS corresponds to a contribution of 
about 66% of a 2 Fa- in the effective a 2 F ab . 

The above results show that, in poly crystalline 
samples and in non-oriented films (where the n band 
dominates the tunneling current) , the inversion of the 
experimental DOS gives mainly information about 
a 2 F n (uj), and thus leads to underestimate the strength 
of the E2 g phonon mode that is essential for supercon- 
ductivity. This might explain the results reported in 
Ref. [12], where the dominant role of the E^g phonon 
mode was questioned. Thus, only separate studies of 
very low-temperature ab-plane and c-axis tunneling 
conductances in high-quality single crystals in clean 
limit might allow a quantitative estimate of the EPI. 
These studies should provide a crucial test for the 
results of the two-band model, even though the fine 
structures in the DOS due to the EPI are very likely 
to be experimentally observable only for tunneling 
current parallel to the ab planes. 
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